Observed Effects of a Changing Step-Edge Density on Thin-Film Growth Dynamics 
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We grew SrTiOa on SrTiOa [001] by pulsed laser deposition, while observing x-ray diffraction at 
the (00|) position. The drop AI in the x-ray intensity following a laser pulse contains information 
about plume-surface interactions. Kinematic theory predicts ^ = — 4<r(l — a), so that =j- depends 
only on the amount of deposited material a. In contrast, we observed experimentally that |^| < 
4<t(1 — a), and that depends on the phase of x-ray growth oscillations. The combined results 
suggest a fast smoothing mechanism that depends on surface step-edge density. 
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Particles of high kinetic energy can interact with a sur- 
face via a variety of exotic, atomistic processes. A goal 
of modern surface science is to describe these interac- 
tions, which occur, for example, in thin-film growth by 
energetic particle beams, such as pulsed laser deposition 
(PLD). Despite the widespread use of PLD to grow high 
quality crystalline films (especially functional complex- 
oxide materials Q, 0, Q), a complete atomic-level de- 
scription of the PLD process remains elusive. In PLD, 
a megawatt laser pulse ablates the surface of a target 
material, creating a highly directed and extremely hot 
and dense "plume" of ionic and neutral species [4J. A 
fundamental understanding of the interaction between 
the plume and the film surface would explain why PLD 
can produce atomically smooth films, rivalling thermal 
(molecular-beam epitaxy, chemical vapor deposition) and 
other energetic (sputter deposition, ion-beam assisted de- 
position) growth techniques. The sensitivity of functional 
properties to film smoothness and crystallinity necessi- 
tates this understanding from a technological standpoint. 
A complete model of PLD film growth may require in- 
clusion of energetic processes not usually considered in 
standard models of growth. 

Conventional growth models typically assume that in- 
cident particles land at random positions on the substrate 
surface. The particles then diffuse until they either (i) 
evaporate, (ii) attach to an existing step edge, or (iii) 
run into other particles and nucleate a new island ||- 
Depending on the substrate temperature and specific en- 
ergy barriers, film growth proceeds in one of several well- 
known modes: 3D, layer-by-layer, or step-flow 0,0,0- 
This treatment neglects the kinetic energy of the incident 
particles (Ek <C 1 eV for "thermal" sources), typically 
10s of eV in the PLD plume 01 ■ Other technologically 
important growth techniques such as sputter deposition 
and ion-beam assisted deposition also utilize particles 
with kinetic energies of 10s to 100s of eV, to enhance 
smooth growth. Previous kinetic Monte Carlo-molecular 
dynamics simulations [fSLllCj. and scanning-tunnelling mi- 



croscopy studies of nucleation density 0, , suggested 
that metal ions with energies near 20 eV can "insert" 
themselves into surface islands if they land near the is- 
land edges, during hyperthermal epitaxy. These inser- 
tion events occur during the initial impact and are com- 
plete well before particles come into thermal equilibrium 
with the substrate (r t / lermQ ; ~ 1 ps). The efficacy of 
this insertion mechanism depends sensitively upon the 
step-edge density and thus upon the surface morphology. 
In this letter, we report time-resolved x-ray scattering 
evidence of a qualitatively similar morphological effect 
during SrTi0 3 (STO) homoepitaxy via PLD. 

Our PLD chamber is an integral component of a fully 
featured x-ray diffractometer that is permanently in- 
stalled in the G3 hutch of the Cornell High Energy Syn- 
chrotron Source (CHESS). By controlling the chamber 
oxygen pressure, Fo 2 , substrate temperature, T, and 
laser repetition rate, /, we can select different growth 
modes, while monitoring film growth in-situ with x- 
rays. We grew STO films on STO [001] substrates in 
two growth regimes: 1) Po 2 — 10 mTorr, T — 540°C, 
/ = 0.03 Hz and 2) P ,_ = 0.01 mTorr, T = 750°C, 
/ = 0.1 Hz. We set the laser energy density at the single 
crystal STO target, and the target-substrate distance, to 
provide an average cr « 0.1 monolayers/pulse. To pre- 
pare a Ti02-terminated surface, we etched the substrates 
in buffered NH 4 F-HF [H Q. We annealed some sub- 
strates at 1000° C prior to growth to obtain a smooth 
surface 0,0. 

A key feature of x-ray scattering is that kinematic scat- 
tering theory accurately describes the scattered intensity 
and can directly test proposed growth models. For the 
simple case of momentum transfer normal to the surface 
of a non-miscut crystal, one has, in homoepitaxy |l7j . 
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Here q is the scattering vector, F (q) is the scattering am- 
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plitude of a single layer, 9 n (t) is the time-dependent cov- 
erage of the n th layer, and d is the layer spacing. The first 
term represents the scattering from an ideally terminated 
single crystal. The second term represents the scattering 
from the deposited film. At the "anti-Bragg" (00^) po- 
sition, x-rays scattered from adjacent layers interfere de- 
structively, providing maximum sensitivity to single step 
height fluctuations. By monitoring the anti-Bragg inten- 
sity, one can monitor the layer coverage in real time and 
directly correlate the growth mode with the experimental 
conditions. This technique has seen wide use in reflection 
high-energy electron diffraction 0,H,EH1(|, and x-ray 
growth studies 

mum. 

Under the conditions of regime 1, we monitored the 
anti-Bragg intensity during STO homoepitaxy [dots in 
Fig. ^a)] . Each oscillation corresponds roughly to the 
completion of one monolayer. Each data point repre- 
sents the scattered intensity integrated over a 15s win- 
dow preceeding a laser pulse, up to, but not including, 
the fast change AI associated with the pulse. The films 
initially grow in a quasi-2D layer- by-layer mode 0,13- 
After the first few layers, the growth takes on a 3D char- 
acter. The 3D growth is equivalent to roughening of the 
growing surface, which gradually decreases the oscillation 
maxima. Additionally, a slight variation in the growth 
rate (~ 0.001 monolayer/pulse) across the 1 cm sample 
causes the oscillation minima to increase during growth, 
as scattering from different regions of the surface becomes 
slightly out-of-phase. Accounting for these two factors, a 
fit to the data accurately describes the evolution of the 
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FIG. 1: (a) Anti-Bragg x-ray intensity oscillations dur- 
ing homoepitaxy of ~ 42 ML of SrTi0 3 ; Po, = 10 mTorr, 
T = 540° C, / = 0.03 Hz. The fit (solid line) to the data 
(dashed line) is generated using Eq. Q. (b) Evolution of 
layer coverages 9 n , determined by the fit. Each line repre- 
sents the growth of one ML. Inset: Surface RMS roughness 
o~rms versus film thickness t during growth (both axes in 
lattice units). Dashed line indicates random growth. 



anti-Bragg intensity over the entire 42-layer deposition 
[solid line in Fig. ^a)]. From the fit, we extracted the 
9 n (t) [Fig.^b)]. Each 8 n (t) initially increases linearly as 
pulses of ablated material arrive at the surface. As antic- 
ipated, new layers nucleate before completion of under- 
lying layers, increasing the number of exposed layers as 
deposition proceeds. From the measured 9 n {i), we have 
calculated the time-evolution of the RMS roughness of 
the surface (inset to Fig.Q. As expected, at late times, 
the roughness grows as urms ~*^>/?= 0.5 [25j |. 

While this kind of measurement and analysis provides 
useful information on surface morphology and growth 
modes, much more information can be obtained by con- 
sidering the pulsed nature of the growth. In the im- 
pulse approximation, deposited material lands at ran- 
dom lattice sites on the (partially) exposed layers, in- 
stantaneously changing the layer coverages by A9 n = 
(6 n — 9 n+ i)(T. Calculating AI using A9 n and Eq. 
leads directly to the elegantly simple result 
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where / is the pre-pulse intensity. Equation J2J is the 
first significant new result of this work. Note that ^y- 
is a function of a only and is independent of both the 
surface morphology and the growth mode. It thus es- 
tablishes a rigorous criterion for determining whether a 
system deviates from simple random deposition during 
pulsed homoepitaxy. 

To illustrate the significance of this result, we have 
simulated the anti-Bragg intensity during pulsed deposi- 
tion of 1.2 monolayers of material onto a perfectly flat 
substrate. Fig. |21 illustrates both the instantaneous in- 
tensity drop AI associated with every pulse and the sub- 
sequent smoothing due to thermal relaxation. The ther- 
mal smoothing is less obvious during the first half of the 
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FIG. 2: Simulated quasi-2D growth of 1.2 ML, at growth 
rate \ = 12 pulses/ML. 
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oscillation, because the overlayer coverage is small, so 
most deposited material lands on the substrate. Thermal 
smoothing becomes visible when the overlayer exceeds 
half-coverage, because the majority of deposited material 
lands on the overlayer. Note that AI is always negative. 
Random deposition always roughens the surface. The in- 
tensity oscillations associated with layer-by-layer growth 
require interlayer transport to smoothcn the film. 

To test the simple prediction of Equation J3J), we mea- 
sured —f- during STO homoepitaxy via PLD, under the 
conditions of growth regime 2, in which films grow in a 
near steady-state layer- by- layer mode 0. To obtain 
the required time resolution, we used a multi-channel 
scaler (MCS) with 1 ms dwell time per channel, and 
10000 channels per pass. The MCS triggered the laser 
at the midpoint of each pass, so that we recorded 5s 
of detector intensity before and after each laser pulse. 
During the 5 s prior to the laser pulse, interlayer diffu- 
sion has essentially concluded, so that the intensity is 
nearly constant. For analysis, we normalized each MCS 
pass by the intensity / integrated over the pre-pulse re- 
gion. Since =}- should be independent of the surface mor- 
phology, we averaged all normalized MCS passes from 
one deposition. We fit a straight line through the pre- 
pulse data, and an exponential through the post-pulse 
data 0|- Evaluating the difference between the fits at 
t = 5s gives = -0.108 ± 0.003. From the growth os- 
cillation data, we find a — 0.078 monolayers/pulse, and 
— 4cr(l— a) = —0.287, clearly indicating a discrepancy be- 
tween theory and experiment (Fig. 13}. STO homoepitaxy 
via PLD is therefore more complex than the pulsed, ran- 
dom deposition implicit in the impulse approximation. 
This is the second new result reported here. 

The discrepancy widens dramatically when we consider 
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the effect of the surface morphology. Over the cycle of 
one intensity oscillation, the step-edge density varies as 
islands nucleate, grow and coalesce [27j. For purposes of 
discussion and analysis, we have assigned an angle <f> to 
specify the phase of the intensity oscillations (Fig. |2J). 
It is reasonable to assume that the density of step edges 
is lowest near intensity maxima (cf> = 0°), increases as 
islands nucleate and grow (0° < cf> < 180°), reaches a 
maximum (cf> sw 180°), and then decreases as the islands 
coalesce (180° ?S 4> < 360°). Therefore, we interpret <j> to 
be the phase of step-edge density oscillations. 

The surface independence of in pulsed homoepi- 
taxy equates with 0-independence. To test this pre- 
diction, we averaged the normalized MCS passes at a 
particular (f>. A visual inspection of the data at three 
phases clearly reveals a morphological dependence to ^ 
[Fig. Ufa)]. While the shapes of the intensity curves agree 
with qualitative predictions of the model [Fig. [5J, the 
magnitude of A£ clearly exhibits a phase dependence, 
with |^| < 4<t(1 — a) at all </>. This is the third new 
result reported here. Fig. 0Jb) illustrates quantitatively 
the dependence of =£■ on <f>. 

The particular morphological dependence exhibited in 
Fig. 01 is consistent with an energetic smoothing mecha- 
nism, the activation of which depends on step-edge den- 
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FIG. 3: Average of all normalized MCS passes for sample 
grown at Po 2 = 0.01 mTorr, T = 750° C, / = 0.1 Hz. Solid 
lines represent linear (t < 5 s) and exponential (t > 5 s) fits. 



FIG. 4: (a) Average of normalized MCS data at three growth 
oscillation phases <j>. (b) Measured at all experimental <j>. 
Dotted lines represent theoretical upper and lower bounds. 
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sity. The data (e.g., Fig. |3J demonstrate that thermal 
smoothing occurs on time-scales on the order of seconds. 
More importantly, they reveal the existence of another 
smoothing mechanism, instantaneous relative to the 1ms 
binning times used here. Interestingly, the magnitude 
of most closely approaches the theoretical value near 
4> Kt 120°, and not at growth oscillation maxima (<fi = 0°), 
where one expects the lowest density of step edges. How- 
ever, the smallest step-edge density can occur at <p > 0° 
if small islands nucleate on top of a rapidly coalescing 
underlayer. The maximum step-edge density would oc- 
cur just prior to coalescence, which could coincide with 
the maximum deviation of ^ (4> ~ 230°). 

In summary, the non-thermal smoothing mechanism is 
very fast, and its efficacy is correlated with the step- 
edge density. Presumably, the physical nature of the 
mechanism is more complex than the insertion mecha- 
nism suggested for noble metal systems. However, step 
edges may provide chemically favorable incorporation 
sites. Highly kinetic plume species landing near these 
sites would rapidly overcome potential barriers to film 
incorporation. An alternative explanation would be a 
<p— dependent deposition rate, a possibility we are cur- 
rently investigating. While this interpretation could ex- 
plain the variation of =p with cf>, the analysis of Fig. 
nonetheless supports rapid smoothing. 

This type of time-resolved x-ray measurement and 
analysis is not restricted to PLD but applies equally 
well to all types of pulsed growth techniques where the 
impulse approximation applies. Thus, any deposition 
beam that can be chopped on time-scales short relative 
to the characteristic time-scale for thermal transport can 
be studied. Examples of potential systems include elec- 
trostatically switched sputter deposition, hyper-thermal 
metal-ion beams, and pulsed supersonic beams. One can 
generalize Eq. 0) to apply this analysis to heteroepitaxy. 
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